Combining multigrid and wavelet ideas to construct more efficient iterative multiscale algorithms
نویسندگان
چکیده
It is shown how various ideas that are well established for the solution of Poisson’s equation using plane wave and multigrid methods can be combined with wavelet concepts. The combination of wavelet concepts and multigrid techniques turns out to be particularly fruitful. We propose a modified multigrid V cycle scheme that is not only much simpler, but also more efficient than the standard V cycle. Whereas in the traditional V cycle the residue is passed to the coarser grid levels, this new scheme does not require the calculation of a residue. Instead it works with copies of the charge density on the different grid levels that were obtained from the underlying charge density on the finest grid by wavelet transformations. This scheme is not limited to the pure wavelet setting, where it is faster than the preconditioned conjugate gradient method, but equally well applicable for finite difference discretizations.
منابع مشابه
Using Natural Wavelet Bases and Multiscale Stochastic Models for Tomographic Reconstruction
We use a multiscale natural pixel type representation of an object, originally developed for incomplete data problems, to construct nearly orthonormal basis functions. The coefficients of expansion of an object in these basis functions are obtained as the 1-D wavelet transform of the (strip integral) projections of the object. This enables us to formulate a multiscale tomographic reconstruction...
متن کاملAdaptive Algebraic Multigrid
Efficient numerical simulation of physical processes is constrained by our ability to solve the resulting linear systems, prompting substantial research into the development of multiscale iterative methods capable of solving these linear systems with an optimal amount of effort. Overcoming the limitations of geometric multigrid methods to simple geometries and differential equations, algebraic ...
متن کاملMultiscale Algorithm for Atmospheric Data Assimilation
We propose a novel multiscale algorithm for the problem of model assimilation of data. The algorithm allows one to efficiently perform optimal statistical interpolation of observed data from a given forecast wf and vector of observations wo. The core of the new approach is a combination of two multiscale tools: a multiresolution iterative process and a multigrid fastsummation technique. Our app...
متن کاملMultiscale solvers and systematic upscaling in computational physics
Multiscale algorithms can overcome the scale-born bottlenecks that plague most computations in physics. These algorithms employ separate processing at each scale of the physical space, combined with interscale iterative interactions, in ways which use finer scales very sparingly. Having been developed first and well known as multigrid solvers for partial differential equations, highly efficient...
متن کاملApproximate cyclic reduction preconditioning
We discuss an iterative method for solving large sparse systems of equations. A hybrid method is introduced which uses ideas both from ILU preconditioning and from multigrid. The resulting preconditioning technique requires the matrix only. A multilevel structure is obtained by using maximal independent sets for graph coarsening. Ideas from [20], [22] are used to construct a sparse Schur comple...
متن کامل